Kondo hole in one-dimensional Kondo insulators 
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Properties of a non-magnetic impurity in Kondo insulators are investigated by considering a 
one-dimensional Kondo lattice model with depletion of a localized spin. The ground state phase 
diagram determined by the Lanczos method shows that the magnetic moment is more stable than 
in ordinary metals. Temperature dependence of impurity susceptibilities is also studied by using 
the finite temperature density-matrix renormalization group. 

PACS numbers: 75.20.Hr, 75.30.Mb, 71.27.+a 



I. INTRODUCTION 

The single-impurity Anderson model fll| was intro- 
duced to study the condition for formation of localized 
moments in metals. The self-consistent Hartree-Fock 
method showed that there is a transition between the 
magnetic state and the non-magnetic state at Uc ^ ^/ Pd, 
where U is on-site Coulomb energy of a virtual bound 
state, and pd is its density of states at the Fermi level. 
But as a result of the Kondo effect [3| , the sharp tran- 
sition disappears and becomes a cross-over. That is, at 
high temperatures higher than the Kondo temperature, 
there is the region where thermodynamic properties are 
well described by the localized magnetic moment. At 
zero temperature, however, the magnetic moment disap- 
pears by the Kondo effect leading to the singlet ground 
state which is obtained by K.Yosida Q- 

In metals, gap-less excitations around the Fermi en- 
ergy play an important role to the Kondo effect. On the 
contrary, in spin-1/2 anti-ferromagnetic (AF) Heisenberg 
two-leg ladders, which have spin gaps and short range sin- 
glet correlations, a non-magnetic impurity, i.e. depletion 
of a spin, induces a local magnetic moment around the va- 
cant site because the impurity breaks a dimer singlet pair 
of the resonating valence bond (RVB) state and leaves an 
unpaired spin _4., Aj- In fact, the magnetic moments ex- 
hibit a Curie law susceptibility, whose Curie constant is 
proportional to the concentration of the impurities [g- 
Correlations between the magnetic moments are stag- 
gered and rapidly decaying with the impurity distance. 
Experimentally, when concentration of the impurities ex- 
ceeds a critical value, the spin gap disappears and it leads 
to an instability of the spin-liquid state towards an AF 
ordered state 7]. On the other hand, theoretically, the 
critical concentration seems to be zero for a bipartite lat- 
tice and thus any finite concentration of impurities leads 
to the ground state with AF quasi long range order and 
gap-less excitations. 

An impurity in the Kondo Lattice (KL) model which 
describes heavy fermion system is introduced as a de- 
pletion of a localized spin similarly to the impurities in 
the two-leg spin ladder. However, since the KL model 



includes both itinerant electrons and localized spins in 
contrast to the spin ladder systems, after the depletion 
one should ask necessary condition for formation of the 
localized moment on the impurity site. The model Hamil- 
tonian for general dimension d is defined as 
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where s„ := '^g g, ^as^s'cl^gCns' and n is d-dimensional 
vector and n^ indicates the impurity site. This model 
represents a single impurity and the surrounding KL 
model. It can be viewed also as the KL model without 
the spin at n — Ud'. the KL model with a Kondo hole. 
When the KL model at half filling, so-called Kondo insu- 
lator, is considered we may expect that the existence of a 
spin gap and a charge gap will lead to similar properties 
as the two-leg spin ladder in certain situations. In fact, 
by substituting non-magnetic Lu-ions into a Kondo in- 
sulator, Ybi_a:Lu2:Bi2, Curie-like susceptibilities are ob- 
served at low temperatures for small concentrations |3]- 
This fact suggests that localized magnetic moments are 
induced by the substitution of the non-magnetic Lu-ions. 
Correspondingly, the optical gap of the Kondo insulator 
YbBi2, is partially filled by the substitution Q- 

Formation of impurity bands in Kondo insulators is 
discussed by Schlottmann and his coworker by using the 
periodic Anderson model by employing the 1/d expansion 
for the self-energy of the host sites for which the second 
order perturbation is used [I^l • Concerning the impurity 
site. Coulomb interaction is neglected, corresponding to 
e — U — in our notation. 

In this paper, we study the condition for the presence 
or absence of localized moments in Kondo insulators and 
properties of these moments. For this purpose, we will re- 
strict ourselves to the one-dimensional (ID) case (Fig.^, 
setting the hopping parameter i„^„' = —t = —1 {i.e., 
measuring all energies in units oft). In the one-dimension 
we can use various powerful methods like the numerical 
exact diagonalization and the density matrix renormal- 
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FIG. 1: The one-dimensional Kondo lattice model with an 
impurity site, Ud- 



ization group. Various properties of the pure ID Kondo 
insulator: 



H^ 



(n,n') ,s 
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have been studied extensively lll|. In particular, mag- 
nitudes of the charge gap and the spin gap are obtained 
numerically [l2 |. The spin-spin correlation functions de- 
cay with oscillations, the RKKY oscillations, and its cor- 
relation length ^ are also determined p^. 



II. PROPERTY OF THE GROUND STATE 

In order to determine the magnetic or non-magnetic 
ground state, the total number of conduction electrons 
iVtot and the total spin S^^^ are studied. Since the sur- 
rounding Kondo insulator has the totally singlet ground 
state when the impurity site is neglected, we expect that 
there is one-to-one correspondence between A'tot and the 
existence or absence of the magnetic moment around the 
impurity site. Let us start from the strong coupling limit, 
J/t ^ 1. To make the following discussion clear we as- 
sume a finite J and t ^ 0. Then the impurity site is 
independent from the rest and it is sufficient to think 
only four states for the impurity site: the unoccupied 
impurity state {S^^^ — 0, A^tot — L ~ 1), the doublet 
(^t^Qt = ±^,Ntot — L), and the doubly occupied state 
(Sj^Qj = 0, A'tot = -^ + 1), where L is the number of the 
lattice sites including the impurity site. In this case, it 
is clear that these states are distinguished with the to- 
tal number of conduction electrons of the ground state 
(Aftot)gs- For a finite t various hopping processes take 
place. However, these processes keep the number of con- 
duction electrons and total spin of the system. Since 
the quantum numbers mentioned above define a disjoint 
subspace, we can investigate the problem of formation 
of magnetic moment just by looking at the total num- 
ber of electrons at finite J/t. The quantum number of 
the ground state is obtained by comparing the lowest 
eigenvalues with the Lanczos method for different sub- 
matrices. We place the impurity site at the center and 
use the open boundary conditions. 

There is an electron- hole symmetry for the present sys- 
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FIG. 2; Phase diagram of {Ntot}gs{U,e) obtained from the 
Lanczos method using the same L = 7 for J = 1.2, 2.4, 3, 10. 
Inset :comparison between the results by the Lanczos method 
(dashed lines) and the fourth order t/J perturbation (solid 
lines) for J = 10. 



S„ -^ — S„, the Hamiltonian (Eq. ^ is reduced to 

{n,n'),s n^Ud 

■ Uc^^^Cnd]C^^iCnal 



-i^ + U)Yci^sCnas 



+2e + U, 



(3) 



where the last term 2e + U can be ignored. Therefore 
the physical properties for a set of parameters (C/, e) are 
identical to the system with ([/, —e — U). In particular, 
since n„s goes to 1 — n„s under the transformation, this 
symmetry leads to 



(Mot)gs(f/,e) = 2L-(Mot)gs(C/, 



U). 



(4) 



tern. By using the transformation, c„s — > (— l)"cjj^ and 



In the strong coupling limit ,J — oo, the energies of 
-^tot — L — \,L,L + 1 are determined as 0, e, 2e -I- U . In 
particular when e = f7 = 0, these three states are de- 
generated. We define such a point as the tricritical point 
[UcT^c) {i-e., {Uc,ec) = (0,0) at J = oo). In general, 
from the electron-hole symmetry one can show that the 
tricritical point {Uc, Cc), if it exists, must be on the sym- 
metric line ec = —Uc/2. 

Figure [3 shows that as decreasing J, the domain of 
the doublet (iVtot)gs(C^, e) = L {i.e., the magnetic state) 
is extended. The transition lines of J = 10 (the solid lines 
in Fig. [SI are very close to those of the strong coupling 
limit. In the strong coupling case, J/t ^ 1, one can 
use the perturbation theory to estimate the boundaries. 
In the inset of Fig. [3 comparison is made between the 
fourth order perturbation results and the numerical ones 
for J = 10. The agreement is satisfactory. 

We comment on the finite-size effect on the phase 
diagram. The numerical results shown in Fig. [21 are 
for a fixed system size, L — 7. When J > 2.4, the 
system-size dependence of {Ntot)gs{U, e) is very small and 




FIG. 3: The finite size effect on tlie pfiase diagram obtained 
from tlie Lanczos metliod for J = 1.2, t = 1 with varying 
L = 3,5,7. 



J 


Uc 


1.2 

2.4 

3 

oo 


-5.96 <Uc< -4.54 

-1.58 ±0.02 

-0.82 ±0.02 





TABLE I: Tricritical point Uc{J) estimated from the Lanczos 
method, tc is given by —Uc/2. 



has reached a convergence within L < 7. But, when 
J becomes smaller the system-size dependence becomes 
more significant. Fortunately, however, the finite-size ef- 
fect shows an alternating behavior as shown in Fig. O 
From the system-size dependence, we may conclude that 
-5.96 < Uc < -4.54 for J = 1.2. Table HI summarizes 
the results of the tricritical points for various J, where 
the errors are estimated from the finite-size effect. 



III. TEMPERATURE DEPENDENCE OF 
SUSCEPTIBILITIES 



In order to study the process of formation of the mag- 
netic moment as a function of temperature, we calculate 
susceptibilities. The values of susceptibilities at zero tem- 
perature depend on the nature of the ground state. At 
low temperatures, susceptibilities diverge like a Curie-law 
if the ground state is magnetic, while converge to zero if 
the ground state is non-magnetic. 



1. definition of susceptibilities 

With a global magnetic field h, we define total and 
impurity susceptibilities in the following way. 



Xtot :- 

Ximp • = 



9{S^, 



dh 
dh 



(5) 
(6) 



where the z-component of the total spin is given by 



C2 



(E 



n=ind 



^n + Z^r, 



The system consists of 



-'tot 

the L — 1 sites of the Kondo insulator and one impurity 
site. Therefore, in order to extract the impurity effect, 
we define another impurity susceptibility: 
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where Xuk i^ defined by Xtot of the ID Kondo insulator 
without any impurity (Eq. |2J). Of course, the contribu- 
tions to this susceptibility from the Kondo insulator parts 
are zero at zero temperature, because both of them have 
the singlet ground states. 

These susceptibilities are calculated by using the 
finite-temperature density-matrix renormalization group 
(finite-T DMRG) jl3| method as follows. From the maxi- 
mum eigenvalue Am of the quantum transfer matrix T/\/ , 
the partition function Zkl of the ID Kondo insulator 
without any impurity (Eq. [21 is given by 
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where M is the Trotter number. The last equality is 
obtained in the thermodynamic limit L — > oo. Following 
S.Rommer |l5| . we obtain the partition function Z of 
the model Hamiltonian (Eq^l by the quantum transfer 
matrix of the impurity site 7^™''. 

Z = tre-'^^= lim trTAf^/^-i^imp 



M^oo 



'AI 



= lim A^f-i(vI/L^|T™P|vI.f,), wheni^oo(9) 

where Xm and the normalized eigenvectors (^]{,/|, |^m) 
are obtained from the finite-T DMRG for a given trot- 
ter number M. The susceptibility is obtained from the 
free energy F or the magnetization m numerically as 
—SF/^Y' or Sm/Sh under a small magnetic field Sh. Typ- 
ical truncation errors in the infinite-system algorithm of 
finite-T DMRG method are 10~^ at the lowest tempera- 
ture kT = 0.12 for the Trotter number M = 200. 

In addition, non-local susceptibilities can be defined 
by applying a small magnetic field on a site, for exam- 
ple at the impurity site, and measure magnetization at 
different sites. We will use such susceptibilities in lat- 
ter discussions to investigate how effects of the impurity 
moment extend to the bulk of the Kondo insulators. 



2. Result of susceptibilities 

In the strong coupling limit, J — oo, it is simple 
to calculate the susceptibihties, because the magnetic 
moment is completely localized at the impurity site: 
(0) _ ^r ^\Q x+le-'f'l^ , ,,(0) .,(0) (0) 



Atot 

/3 



(L - l)/3 



+ Xi: 



Ai] 



Xd 



e3J;3/4_|_4^g-j/3/3 ' A-imp ' A.imp 

These susceptibilities at low temper- 



'2(2+e'3^+e-/3(=+C')) 

atures show a sharp transition at the boundaries around 
the magnetic solution. In the magnetic region, the ra- 
tio of the susceptibilities to the Curie law susceptibility 
(xc = l/4fcT) converge to unity in the limit of T — > 0: 

(0) / (0) ; T 

At high temperatures T ^ J and f , susceptibilities 
of any J/i show good agreement with the susceptibili- 
ties of the strong coupling limit, since correlations be- 
tween different sites are negligible. In fact, the high- 
temperature expansion of the susceptibility %, which is 
obtained from the second order perturbation of i/J, is 

given by Ximp - x|mp = -i^/3V4 • 

On the other hand, at low temperatures the values of 
susceptibilities depend on the ground state as we men- 
tioned at the beginning of this section. Figure ^ shows 
Xd/xc for various J and ({/, e) on the symmetric line: 
t/ = 8,0, -8 and e == -C//2. First we note that Xd/xc 
for U — Q and J = oo is special and temperature inde- 
pendent. This is due to the fact that t/ = and e = is 
the tricritical point in the strong coupling limit. Except 
for f/ = and e = case, all lines are pointing towards 
either zero or unity in the zero temperature limit. These 
behaviors are consistent with the values of the tricrit- 
ical point in Table Q] estimated from the ground state. 
However the temperature dependence of Xd/xc is quite 
different for C/ = 8 and 0. 

To understand the difference, we plot Ximp/xc in 
Fig. |S1 One can see that Ximp/xc approaches to a finite 
value for C/ = 8 and while it goes to zero for C/ = — 8. 
Furthermore the finite value seems to depend on U and J. 
As can be seen from the definition of ximp (Eq. El , only 
the induced moment at the impurity site is measured for 
this susceptibility, while the entire moment induced in 
the system is measured for Xd (Eq. 0). Therefore the 
different limit Ximp/xc may be understood as due to the 
difference in the spatial spread of the induced magnetic 
moment. Naturally, we expect that the induced moment 
has a larger amplitude at the impurity site for larger U 
{> Uc) and a wider spread for smaller J since the corre- 
lation length of the Kondo insulator gets longer as J is 
decreased [l^. We may conclude that the limiting value 
of Ximp/xc is larger for larger U (> Uc) and larger J. 
This is exactly what we observe in Fig. [S] 

In order to investigate spatial spread of the magnetic 
moment, we calculate the susceptibilities 



Xs 



5hd 



(10) 



where hd is the local magnetic field applied on the impu- 
rity site only. 
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FIG. 4: Ratio of the susceptibility Xd to the Curie suscep- 
tibility l/4A:r calculated by the finite-T DMRG method for 
J = 3, 2.4 and for e = -C//2; [/ = 8, 0, -8. 




FIG. 5: Ratio of the susceptibility Ximp to the Curie suscep- 
tibility l/4fcr calculated by the finite-T DMRG method for 
J = 3, 2.4, 1.2 and for e = -C//2; [/ = 8, 0, -8. 



When a magnetic field is applied on one site in the 
Kondo lattice model, the induced magnetic moments 
show the RKKY oscillatory behavior 

(s„) ex (-l)("-"'')exp(-|n-,^o|/C), (H) 

where ^ is the correlation length. Thus we expect that 
the non-local susceptibilities Xs„ also show a similar be- 
havior. 

It is interesting to estimate the correlation length from 
the ratio between Xs^ +i and Xs^ through the relation 



= -log(-Xs„ ,i/Xs„J. 



(12) 



Figure ini shows that values of ^ at low temperatures esti- 
mated from the above relation are close to the correlation 
length of the ID Kondo insulator evaluated at T = by 
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. 6: The correlation length ^ estimated from Xs,i and 
^1 for J — 1.2, 2.4, 3. The correlation lengths of the bulk 

Kondo insulator obtained by the zero-temperature DMRG 

method [13] are also shown. 

the zero-temperature DMRG ^^. The value of ^ in- 
creases with decreasing J, which is consistent with the 
divergence of the correlation length at J = 13]. 

Finally we would like to comment on the non- 
monotonic temperature dependence observed in Fig. 4. 
The decrease of the ratio of the susceptibility to the Curie 
susceptibility around T = 2 with lowering temperature is 
due to formation of conduction band. This is a general 
trend which is shown by the above second order pertur- 
bation theory with respect to t/J. On the other hand, 
the upturn around T = 0.5 ~ 0.8 may be traced back to 
the effect of formation of the Kondo singlet states on both 
sides of the impurity. In fact, this idea is consistent with 
the fact that the minimum occurs around the spin gap 
temperature of the Kondo insulators, however, further 
investigation is necessary to draw a definitive conclusion. 



IV. CONCLUSIONS 



In contrast to the naive expectation, the region of sta- 
ble magnetic moment in the (U, e) plane is wider for 
smaller J. It is surprising that the region of stable mag- 
netic moment extends to the region of negative U for a 
finite J. The key to understand this result is the exten- 
sion of the formed magnetic moment. Longer correlation 
length of the Kondo insulator with smaller J means larger 
spread of the magnetic moment. It means that the ef- 
fect of Coulomb repulsion (attraction for negative U) is 
weakened for smaller J. In the weak coupling case, the 
correlation length diverges exponentially as a function of 
J. On the other hand, the charge gap which is linear in J 
[l^ acts in the Kondo lattice parts as a positive effective 
potential. 

It is interesting to compare the effects of depletion of 
localized spins between the two-leg spin ladder and the 
one dimensional Kondo insulator. The both system have 
spin gaps as a common feature. For the spin ladder, 
depletion of a spin always leads to formation of a localized 
moment around the depleted site. For the Kondo chain, 
formation of a magnetic moment depends on ([/, e) . 

A finite concentration of depletion may lead to quasi 
long range order characterized by gap-less spin excita- 
tions in the spin ladder. On the other hand for the 
Kondo insulator, one can speculate that there may be 
two cases. In one case, the ground state may be non- 
magnetic, while in the other case there may be quasi 
long range anti-ferromagnetic ordering. This will be an 
interesting future problem. 

Finally let us consider the effect of non-magnetic im- 
purities in real three dimensional Kondo insulators. We 
believe that the first conclusion that the formation of 
magnetic moment is generally easier in Kondo insulators 
than in ordinary metals is also valid in higher dimensions. 
It is also probable to consider that there may be the two 
cases, non-magnetic and magnetic ground states for finite 
concentrations of non-magnetic impurities in the Kondo 
insulators in realistic three dimension. However, estima- 
tion of the critical concentration for each case needs more 
elaborate investigation and clearly beyond the range of 
the present paper. 



In this paper we have studied the problem of formation 
of magnetic moment when a Kondo insulator is doped 
with an impurity without the localized spin. One dimen- 
sional Kondo lattice or Kondo chain is treated in this 
article. The ground state phase diagram is obtained by 
the Lanczos method. Temperature dependence of the im- 
purity susceptibilities calculated by the finite-T DMRG 
are consistent with the phase diagram. 
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